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ABSTRACT 


A Taylor-Couette cell for the investigation of geometrical 
phase in acoustics has been constructed. The inner and outer 
cylinders are made of acrylic. The cell has an inner cylinder 


Pacaus of 9.477 cm, radius ratio of 0.902, and ‘cavity aspect 


ratio of 49.2. Two high performance d.c. motors can rotate 
the cylinders independently. The angular speed of the 
cylinders irS-SsMenteerea » iby .«a ‘*pHeOLO-ineerrupLter. A 


piezoelectric polymer (PVDF) transducer mounted on the inner 
cylinder drives the acoustic field in the cavity anda 6.02 mm 
diameter electret microphone embedded in the wall of the inner 
cylinder acts aS a receiver. Static measurements of the 
acoustic modes in the annular cavity show good agreement with 
theory, though measurements at different angular locations of 
the transducer reveal 2% nonuniformities. The effects of 
nonuniformities could be overcome by rotating the cylinders at 
high speeds. The performance under these preliminary tests 
indicate that an apparatus of this design is suitable for 


investigations of acoustics in rotating flows. 
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I. INTRODUCTION 


The motion of a fluid flow between two concentric, long 
cylinders undergoing differential rotation has been one of the 
most studied problems in hydrodynamics (Donnelly, 1991), with 
however, numerous unanswered questions. Such fluid motion is 
called Taylor-Couette flow. 

The steady state axis-symmetric flow is an exact solution 
of the Navier-Stokes equation in this configuration. This 
flow is simply described as the superposition of solid body 
rotation plus the flow of a vortex line along the common axis 
of the two cylinders. When the inner cylinder is turned at a 
higher rate, such that the inner layers of the fluid are 
moving more rapidly than the outer ones, the centrifugal force 
pulls these inner layers outward which then break into cells 
and circulate. The fluid breaks into horizontal bands. Thus, 
the axis-symmetric flow is, in general, not stable. As the 
angular velocity of the inner cylinder is increased, the 
number of bands increases at first, then they become wavey. 
The velocity of rotation of this wavey pattern approaches a 
numerical constant (which is determined by the geometrical 
parameters of the system) times the angular velocity of the 
cylinder. Finally, for higher values of the angular velocity, 
Or in the configuration in which the outer cylinder is 


rotating in the opposite direction, the flow becomes turbulent. 


The complexity of behavior in this simple system is in 
principle contained in the equations of a viscous fluid. 
However, there does not exist, up to the present, a 
satisfactory method to analyze and explain the complexities 
hidden in the Navier-Stokes equation: instabilities, the onset 
of turbulent flows, and intermittency. 

Consider however, the flow not aS a dynamical system, but 
as a collection of parameters of an acoustic field. In the 
geometrical acoustic approximation, the acoustic field is 
described by an equation analogous to Schrdédinger’s equation 
of a charged particle, with the flow velocity playing the role 
of the magnetic vector potential. Then, it is possible to 
determine global properties of the flow by measurements of the 
acoustic phase. Modal resonance splitting and shifting are 
manifestations of global changes. 

This thesis describes the design, construction, and 
preliminary testing of a Taylor-Couette cell. The flow 
apparatus is designed with maximum flexibility in mind. 
Preliminary testing and experimentation using the cell 
determined its Suitability for future research. 

The organization of this thesis is as follows. Chapter II 
presents the theory for the mode structure of the Taylor- 
Couette cell in its static and differential rotation 
configurations. Chapter III describes the design and 


cONSETUGCEION: CLE hemCceIae Chapter IV describes’ the 


preliminary tests and their results. Chapter V presents a 


Summary and conclusions of this work. 


fy THEORY 


A. EXPECTED RESONANCE FREQUENCIES 

In order to observe the modal frequency split and shift in 
a rotating circular Couette system, we require knowledge of 
the modal structure. The following theory derives those 
frequencies. 

Consider the annular cavity formed by two concentric 
cylinders of inner and outer radii, R, and R,, respectively, 
with top and bottom boundaries located at z=0 and z=L,. If 
all surfaces of the cavity are perfectly rigid, so that the 


particle velocity normal to the surface is zero, then for 


linear acoustics f+ Vp=o at all boundaries, where Vp is the 


gradient of the acoustic pressure. Thus 
(2p a 


ar ) ne : (32 


- aT cag 
(1) 


For a system with circular symmetry, it 1s best to recast 
the linear inviscid wave equation in polar coordinates 


Gp ~ lep (cp) soe cee (2) 


Or? sx Or =r? 9B2.—s z2—se2.-:« BEEZ’ 














where c is the speed of sound of the medium. The technique of 
separation of variables is use to solve Eq. (2). The acoustic 
pressure 1S assumed to be a product of four terms, each a 


function of only one variable, i.e. 


pir ,0ue, 6) — Ri rT) Oi) 7(4)e22"; (3) 


Suosctliue len Of EQ. —(3) 1nto EG. (Z)*yi1elds 








dR , @ZdR , RZd°@ , ROd’Z , p2pgz- o. 


@Z Pal Eohac hd 
dr? gota r°de2 dz? 


(4) 
where k=@/c. Multiplying each term in this equation by 1/RO@Z 
and grouping terms containing z on one side results in 


d?R poe. ey, 1d*e ; -1d?Z 


eo k = a 
Rdr? rRdr r*8de2 Zdz? 








(5) 


tieylefti-hand “side of this equation; “a function of r and G, 
cannot be equal to the right-hand side, a function of z alone, 
unless both functions equal some constant. If we let this 


constant be k,’, the right-hand side becomes 





2 
( + K2}z = 0. (6) 


Repeated performance of the separation of variables technique 


applied to Eq. (4) results in the final set of equations: 














2 
HS + ile = 0, (7a) 
Zz 
Zo nl = 0, (7b) 
Ss peelicias: (is = aie = 0, (7c) 
dia) ene ae 


where m? is an arbitrary separation constant and k,’=k’-k,?. 


The complete general solution to Eq. (7a) is 


8 
Z(z) = A,cos(k,z) + A,sin(k,z), sor 
where A, and A, are arbitrary constants. The solution to Eq. 


(7b) 1s 


MC (SCO) = ON oo, 


where B iS an arbitrary constant and @ iS an arbitrary phase 
angle dependent on the location of the source. Eq. (7¢) age 
known as Bessel’s equation. Solutions to this equation are 
the transcendental functions called Bessel functions of the 


first kind J,(k,r) and second kind Y,(k,r) of order m, 


RUA) = Cid. een) + Oe Cara (10) 


where C, and C, are arbitrary constants. Bessel functions are 


oscillating functions whose amplitudes decrease as their 


arguments increase. The second kind functions Y,(k,r), also 
known as Neumann functions, become unbound in the limit as k,r 
aperoaches Zero. 

Applying the appropriate boundary conditions to Eq. (8) 


results in 
Z(z) = A,cos(k,z), 


where 
(11) 
[oats as eI Ome On oe 
If the acoustic pressure p is to be a single-valued function 
Sesposcitlon, then p(r,0,2,t) must equal "p(n, 04+2n,z;t), which 
restricts the separation constant m to integral values; 
MeOp tlh, 2, «« 


When k,=0 and m=0, Eq. (10) implies the solution 


R(r) = C, = constant. (22) 


When k, 1S non-zero, the boundary conditions applied to Eq. 


(10) results in the equality 


Jae.) 7 Tn (k,, R,) 
¥(k,R,) Yi(K,R,) 


(13) 
where the prime denotes the first derivative with respect to 
kr evaluated at the specified argument. If we let the ratio 


R,/R,= g, then Eq. (13) can be written as 


Tm(KnR,) — In(GkKpR,) 
Y,(k,,R,) Y, (gk,, R,) 


=" (oe (14) 


The roots of Eq. (14), values of k,R, can be solved for using 
a computer for various values of m (see Appendix). If the 
values of the m,, order Bessel functions which cause Eq. (14) 
to equal zero are designated by j,,, then the allowed values 
of k, assume discrete values given by k,, = jJn,/R, - 
Appropriate solutions to Eq. (3) for the case k,, not equal 


to zero yield the normal modes 


VY. (omer 
ya etalon” ers) 


Delt, 0, Ze). See ae en) ; 
Oi ew (15) 


x cos (10 +o) cos (k,z) eIe tant 


where 
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elie Os eee 





ce 
k= = @=0,1,2 
ioe 
a 
| a — et m=0,1,2 
m7 R, 
W="C, 1, 20k 
=n+e#odo., 


(16) 


The allowed frequencies are thus quantized by the relationship 


fan = —VKnn + ke (17) 
Each eigenfunction given by Eq. Gls5)) has its own 
characteristic eigenfrequency, denoted by Eq. (17), and can 
be specified by the ordered integers (@,m,n). The value of 0 
determines the number of nodal surfaces in pressure 
perpendicular to the z axis. The value of m specifies the 
number of radial pressure nodes. The value of n gives the 
number of nodal surfaces of constant r (cylinders) surrounding 
the z axis and concentric with it. 

Solutions to Eq. (3) when k,, equals zero yield the normal 


modes 


Pi,0,0 (r,6, Zy t) = Pig gCOs (KsZ) eset 


where 


C 
(18) 
k,= = ¢=0,1,2, 
Zz 
The allowed frequencies are 
fr = eS 19 
,0,0 QT z ( ) 


corresponding to pure z modes. 


B. MODE STRUCTURE FOR DIESERENTTIAL ROTATION 

Axis-symmetric flows in a Taylor-Couette cell are derived 
for the annular cavity between two infinitely long cylinders. 
However, the expected resonance frequencies derived assumed 
finite length cylinders. A justification for assuming such 
velocity field has been given previously by Andereck et al., 
1983. These researchers noted that if the aspect ratio I, the 
ratio of the fluid height to the gap between the cylinders, at 
the region of measurement iS properly chosen, then the 
boundary conditions existing at the top and bottom of the 
cylinders exerted no observable difference in the flow state 
properties. In other words, boundary conditions at the top 
and bottom of the cylinders can be neglected if measurements 
are made far removed from these boundaries. Many of the 
measurements cited in the above re were conducted with an 
aspect ratio of 30. However, other measurements were made 
with 20SI°$47. Our initial testing is conducted with T=24.6. 
Therefore, we can assume that our region of interest has an 
azimuthally symmetric velocity field. 

Consider waves with high azimuthal numbers m. For these 
waves the velocity field induced by the cylinders’ 
differential rotation varies over distances large compared 
with the wavelength. Within this approximation, the local 


dispersion relation in the moving medium with velocity flow 


u(r) is given by 


10 


Oe Oe ane tr): (20) 


where @(k) is the dispersion relation of a wave with wave 


vector kK in the stationary medium. A particular realization 


of the flow velocity is given by the stable azimuthal 
configuration with a stationary inner cylinder. For flow 


velocities lulec, Eq. (20) can be written to lowest order as 





wo = O(|K+ AU |, 2). (21) 
The correspondence between @ and the Hamiltonian establishes 
the analogy between the magnetic vector potential A and the 


flow velocity, provided 





Z x. ~ka 22 
h a ise 
(Berry et al., 1980; Larraza, 1989) Making the natural 


identification between the wave vector k and the operator 


-iV, Eq. (21) is analogous to Schrédinger’s equation for a 


charged particle in a vector potential A. The geometry of 


Taylor-Couette flow for the azimuthally symmetric case yields 


the eigenfrequency spectrum 


hak 


|mlQR, (23) 








OS alee 
CG R 


oO 
(Larraza, 1989), valid for rotation of the outer cylinder 
only. If the rotational speed Q=0, the states m#0 are doubly 
degenerate. For a non-zero value of the rotational speed the 
degeneracy is broken. Physically, this means that the waves 
travelling in the same direction as the average flow have an 
additional contribution to the sound speed, whereas waves 
travelling in opposite direction to the flow have an effective 
sound speed which is less than the static value. According to 
Eq. (23), there is an asymmetric splitting of the degenerate 
states about the zero-angular speed eigenfrequency. To see 


this, notice that Eq. (23) is the first order expansion of 


mie 
_ | Fo (24) 
SE QR, 
1+ 





|m|QR, 
Cc 


Inoue == with n an integer, degeneracy occurs. For n=l, 


the states m=1 and m=-2 are degenerate with frequency 


ssic 
2R. 





, while the state m=-1 is not degenerate. Physically, 


these states are reached by increasing the rotational speed 2 


and correspond to high azimuthal numbers m. 
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AO Gh Ge TAYLOR-COUVETTE FLOW APPARATUS 


A. MECHANICAL DESIGN 
With flexibility in mind, the cell is designed to allow 
the inner and outer cylinders to be rotated independently. 
Also, rigid boundaries at the top and bottom of the cylinders 
allow for three types of boundary motion: stationary, rotation 
at the angular speed of the inner cylinder, and rotation at 
the angular speed of the outer cylinder. The system is 
designed to be vertically mounted to eliminate flexural strain 
resulting from the weight of the components if mounted in a 
horizontal configuration. 
1. Cylinders 
The inner and outer cylinders are made of acrylic. 
The acrylic is fabricated using extruded, vice cast, methods, 
which introduced nonuniformity effects which will be described 
later. The inner cylinder radius,R, is 9.477 cm and the 
outer cylinder radius, R,, is 10.510 cm. The resulting radius 
ratio is 0.902, which is near the value used in many previous 
Studies of Taylor-Couette flow. (DiPrima and Swinney, 1981) 
The solid boundaries used to enclose the annulus create a 
resonant cavity height, L,, of 50.838 cm, resulting in a 


cavity aspect ratio of 49.2. 
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2. Support Structure and Bearings 

With the exception of the cylinders, drive motors and 
their associated drive belts, the entire assembly, including 
Support structures, are made of aluminum. The inner cylinder 
is supported by a length of solid round stock, with a diameter 
of 5.08 cm, which also provides axial alignment for the 
bearings. 

Radial bearings are used to support rotating 
components since radial rigidity is of paramount importance. 
The type of radial bearings chosen are able to handle compound 
loads. Thus, while providing rigidity, they are also able to 
support load resulting from the weight of the assembly. 

3. Upper and Lower Boundaries 

Aluminum rings provide an acoustically rigid boundary 
COmcat aon. The rings can either be stationary, or rotated 
with the same angular speed as the cylinder to which they are 
attached. The rings are machined to fit flush against the 
cylinder they are connected to while providing as narrow a gap 
aS possible near the other cylinder. While a tight cavity is 
desirable to set up standing waves, some gap between the rings 
and the cylinder walls is necessary to eliminate frictional 
rubbing which results in torque disturbance and noise pick-up 
on the receiver. The average gap between ring and cylinder 


wall -ism0-0d= ine 
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4. Motors and Drive Belts 
Two high performance, fractional horsepower, disc- 
armature, dc motors are used to rotate the cylinders. These 
motors, PMI U-series model #U12M4H, are designed for excellent 
speed control and constant torque delivery. The motors are 
connected to the geared cylinders via a stainless steel cable 
reinforced, polyurethane, toothed drive belt. 
5. Miscellaneous 
The entire assembly iS mounted on a 2 cm thick 
aluminum baseplate meaSuring 50 cm by 62 cm. Two vertical 
aluminum cylinders Support a horizontal crossbar which 
provides axial alignment of the components at the top of the 
assembly. Figures 1, 2, and 3 show the construction of the 
apparatus as well as close-ups of the drive assembly and top 


cylinder mounting. 
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Taylor-Couette cell flow apparatus 


Figure 1 
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B. ELECTRICAL DESIGN 
1. Transduction 

Several Constraints; mechanical as well as 
theoretical, led us to choose a piezoelectric polymer film as 
micur Source —OF “.Cdvity ) “excitation. The polymer is 
polyvinylidene fluoride (PVDF). PVDF is sold commercially in 
sheet form, resembling a piece of thick plastic wrap. For our 
application, a 40Um thick strip was cut 2.0 cm wide and 50.8 
emelong tO match the length of the Cavity. For a-constant 
driving voltage, PVDF acts as a constant displacement source. 

Our ultimate goal is to investigate changes in 
azimuthal mode structure due to the rotating flow. 
Consequently, it is not necessary to excite radial and 
vertical modes in order to observe the frequency split or 
shift. Choice of excitation frequency and source placement 
can limit excitation of these modes. 

Since the PVDF strip 1S a constant displacement 
source, its width has to be less than one-half wavelength for 
the range of frequencies of interest. The upper frequency 
limit for preliminary testing was chosen to be 8kHz. Similar 
reasoning as that used for strip width was also applied to 
Strip placement. With the strip extending across the length 
of the cavity, those vertical modes requiring a node at the 


mid-height of the cavity would be weakly excited, if at all. 


B. 


2. Receiver 

A 6.02 mm diameter electret microphone is embedded in 
the wall of the inner cylinder so that its face is flush with 
the surface of the cylinder. A silicone gel is placed around 
the microphone to provide a tight seal. The desire to 
minimize the number of penetrations through the cylinders 
during preliminary testing led us to use only one microphone 
placed at the vertical center of the cavity, 25.4 cm from 
either vertical boundary. However, the use of only one 
microphone and its placement necessarily limits’ the 
interpretation of the results of the preliminary testing as 
will be discussed in Chapter IV. 

3. Angular Speed 

The design uses a photo-interrupter to measure the 
angular speed of the outer cylinder. A small tab, mounted at 
the top of the outer cylinder, gates the output of the photo- 
interrupter once each revolution. The interrupter signal 
output is fed to an HP 5384A frequency counter which counts 
the number of events occurring within a .1l, 1, or 10 second 
window. Digital readout could either be in hertz or the 


average period of revolution. 
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HA RESULTS OF PRELIMINARY TESTS 


To assess the capabilities and limitations of our initial 
design, we performed two preliminary tests. The purpose of 
these tests was to verify our theoretical analysis of the cell 
and to determine the magnitude of errors caused by 
nonuniformities. 

With the inner and outer cylinders fixed, the output of an 
HP 3562A dynamic signal analyzer was used to excite the 
cavity. A 5 Vrms periodic chirp source output, a fast sine 
sweep over the selected frequency span that repeats with the 
Same period as the time record, was used. The output of the 
embedded electret microphone was fed back to the signal 
analyzer and its power spectrum obtained. 

Frequencies corresponding to clearly discernable peaks in 
the power spectrum were recorded. The power spectrum was 
analyzed from d.c. to approximately 3.6 kHz in 400 Hz windows, 
Mmavang 0.5 Hz resolution. Within this frequency span, 
approximately sixty theoretical eigenfrequencies are possible. 
To aid in the resonant frequency cataloguing, an HP 3314A 
function generator was used in conjunction with the dynamic 
Signal generator to manually sweep through a resonance and 
observe the phase change on an oscilloscope. 

The results of these measurements are shown in Figure 4, 


which is a graph of observed resonance frequencies versus the 


Zi 


corresponding theoretical frequencies. The solid line has a 
slope of one. Ideally, the data points should fall on this 
line. The degree to which they do not is shown in Figure 6, 
which shows the magnitude of percent error. Generally, the 
results are quite good. The agreement is typically within 1%. 

A number of modes were easily observable indicating the 
gap between the boundary rings and the cylinder walls did not 
kill the mode structure. The mode structures of interest, 
those with purely azimuthal dependence, were often masked by 
much stronger, closely spaced vertical modes. Vertical mode 
Structures were not eliminated to the degree we had 
anticipated. This maybe due in part to the way the PVDF is 
attached. The strip is actually composed of two shorter 
Strips, joined at the midpoint of the cavity. The joint 
construction may result ina smaller displacement of the PVDF 
at the joint as compared to the rest of the strip, therefore, 
allowing vertical modes to be excited. Further work in this 
area is needed. 

The use of one microphone, located at the center of the 
cavity, is not sufficient to distiguish between closely spaced 
modes. A moveable probe microphone might be better suited for 
determining mode structure. 

Lastly, we did not observe any purely azimuthal modes and 
only about half of the modes represented in Figure 4 have any 
azimuthal component. This again points to the need for 


continued work on the method of transduction. 
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Figure 4 Observable resonance peaks from dc to 
SEC <lglvs 


Abs(Percent Error) 





0 1000 2000 300 4000 
Theoretical Frequency (Hz) 


Figure 5 Agreement between observed resonance and 
theoretical resonance 
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While making the measurements discussed above, we noticed 
that the resonance frequency of a mode depended on the 
relative angular orientation of the inner and outer cylinders. 
We suspect that this is due to nonuniformities in the cylinder 
walls, ring boundaries, and alignment. The magnitude of this 
effect is illustrated in Figure 6. The data show the result 
of measurements of the resonance frequency at six equally 
spaced relative orientations. In Figure 6, the percentage 


variation of the resonance frequencies is plotted as a 


function of the mean resonance frequency. In all cases, the 
percentage variation is less than approximately 1.5%. This 


could account for much of the error shown in Figure 5. 
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Figure 6 Variation of resonance through 360 degrees of 
rotation 
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Vv. SUMMARY AND CONCLUSIONS 


This thesis reports the design and construction of a 
Taylor-Couette cell. The apparatus uses acoustic transduction 
as a means to characterize global properties of the flow. 
Preliminary tests show that the static mode structure differs 
slightly from theoretical predictions, mainly because 
nonuniformities in the cavity. The nonuniformities were 
qualitatively observed by measuring the resonant frequencies 
for different angular positions of the transducer. These 
preliminary tests provide qualitative limitations imposed by 
the apparatus for the applicability of the theory in Chapter 
et. Because of nonuniformities, certain degeneracies are 
asymmetrically removed, mocking the effect of an arbitrary 
flow. Thus, quantitative knowledge of the effect due to 
nonuniformities is essential in order to determine the minimum 
differential rotation for which the effects due to the flow 
become important. 

The acoustic "visualization" reported in this thesis 
constitutes a novel concept, but further development of the 
technique is required. As opposed to standard visualization 
methods, in principle acoustic transduction does not limit the 
Reynold’s numbers for which the technique is useful. But most 
important, because the acoustic phase can provide global 


properties of an arbitrary flow, the investigations proposed 
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in this thesis may lay the foundations for a first-principles 
phenomenological theory of turbulence that relies on a 


complete set of global variables for its description. 


Za 


APPENDIX 


alee 





28 


MODES OF A CYLINDRICAL ANNULUS 
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